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Abstract 

We determine the exact complex-temperature (CT) phase diagram of the ID Z 6 clock model. 
This is of interest because it is the first exactly solved system with a CT phase boundary 
exhibiting a finite- K intersection point where an odd number of curves (namely, three) meet, 
and yields a deeper insight into this phenomenon. Such intersection points occur in the 3D 
spin 1/2 Ising model and appear to occur in the 2D spin 1 Ising model. Further, extending our 
earlier work on the higher-spin Ising model, we point out an intriguing connection between 
the CT phase diagrams for the ID and 2D Z§ clock models. 
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We report here a determination of the exact complex-temperature phase diagram for 
the ID Zq clock model and show how the results give a deeper insight into certain features 
of higher- dimensional models. The idea of generalizing a variable, on which the free energy 
depends, from real physical values to complex values was pioneered by Yang and Lee [|l[] , who 
carried this out for the external magnetic field and proved a celebrated circle theorem on the 
complex-field zeros of the Ising model partition function. The generalization of temperature 
to complex values, performed first for the Ising model ||, has also been quite fruitful, since 
it enables one to understand better the behavior of various thermodynamic quantities as 
analytic functions of complex temperature (CT) and to see how physical phases of a given 
model generalize to regions in CT variables. Indeed, a knowledge of the complex-temperature 
phase diagram of a model for which the free energy is not known provides further constraints 
to guide progress toward an exact solution. A classic example relevant here is the 3D Ising 
model; the main purpose of the present work is to present some exact results for a ID 
model which elucidate an important feature of the complex-temperature phase diagram of 
the 3D Ising model. Some early works on CT singularities include Refs. [|]-|| studying 
partition function zeros, and Ref. 0, motivated by the effect of these singularities on 
low-temperature series. The continuous locus of points in the complex-temperature plane 
where the free energy is non-analytic is denoted B. For the spin models of interest here, 
with isotropic, nearest-neighbor spin-spin couplings, this is a 1-dimensional curve (including 
possible line segments). Parts of B form boundaries of various regions, some of which are 
complex-temperature extensions of physical phases and some of which may have no overlap 
with any physical phase. B may also include arcs or line segments which protrude into, and 
terminate in, certain CT phases. 

In general, phase boundaries of physical systems include intersection points where differ- 
ent parts of the boundaries meet. Examples include (i) the triple point in the (T,p) phase 
diagram for a substance like argon, water, etc. where gas, liquid, and solid phases all coexist, 
and (ii) the two triple points forming the ends of the A line in the phase diagram for helium, 
where, at the lower end of this line, the gas, normal fluid and superfluid coexist, and, at the 
upper end, the solid coexists with the normal fluid and superfluid. Similarly, the complex- 
temperature phase boundary B of a spin system may include intersection points at which 
different curves contained in B meet. We denote the number of curves meeting at such an 
intersection point as n c . Often these can be grouped into pairs, such that two curves meet 
with equal tangents at the intersection point; these are then regarded as a single branch of 
a curve passing through this point. In the terminology of algebraic geometry ||, a singular 
point of an algebraic curve is a multiple (intersection) point of index m if m branches of 
the curve pass through this point. This is thus an intersection point with n c = 2m curves 
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meeting, in m pairs with nt = m equal tangents. For the 2D Ising model on the (homopolyg- 
onal) square, triangular, and honeycomb lattices, and also on the heteropolygonal 3 • 6 ■ 3 • 6 
(kagome) and 3-12-12 lattices, the CT phase boundaries B are algebraic curves and their 
intersection points always have n c = 4 and n t = 2 so that they are regular multiple points 
of index m = 2 (see Tables 2,3 in Ref. ||). This is also true of the Ising model in a nonzero 
external field satisfying the Lee- Yang condition (3H = iix/2\ the CT phase boundaries B on 
the square, triangular, honeycomb, and 3-12-12 lattices all exhibit intersection points with 
n c = 4, n t = 2 which are thus again regular multiple points of index m = 2 |[10|| . In all of 
these cases, the two branches of the curves cross with an angle 8 C = 7r/2, i.e., orthogonally. 

However, in other models, it appears that the respective phase boundaries B have inter- 
sections at finite complex-temperature points with the odd n c = 3. In contrast to the 2D 
Ising model results noted above, which are extracted (in Refs. p], ||, || [II], |12| ) from exact 
solutions (|14| for square lattice; reviewed for general lattices in Ref. JL5[]), up to the present 
time, to our knowledge, none of the models which appear to have finite- K intersection points 
with odd n c has been exactly solved (here, K = J/(&bT)). Instead, the occurrence of such 
points has been inferred from inspection of CT zeros of the partition function calculated on 
finite lattices. For example, for the (spin 1/2) Ising model on the simple cubic lattice, such 
zeros [|], suggest an intersection point at Ui/ 2 — —0.5, where u s = e~ K ^ s2 . At this point, 
a component of B crossing the negative real U1/2 axis vertically meets a line segment lying on 
this axis in a T intersection, so that n c = 3, n t = 2. Since the lattice is bipartite, there is a 
similar intersection point at the inverse position, U1/2 — —2. Similarly, in our studies of the 
spin 1 Ising model on the square lattice we have found, for the largest lattice sizes, indications 
of n c = 3 intersection points at u\ ~ 0.1 + 0.8z, its inverse, and their complex-conjugates 
j~8| . To gain further insight into this feature of complex-temperature phase boundaries, one 



is motivated to search for an exactly solvable model which exhibits this behavior. We have 
succeeded in this; the ID Z 6 clock model provides to our knowledge the first exactly solvable 
model with a finite- K intersection point having odd n c , and we analyze it here. 

The Z N clock model is defined at temperature T by the partition function Z = J2{e n } e~^ n 
with f3 = (ksT)" 1 and 

H =Y,E(-£ n + t n ,) (1) 

(nn') 

where the interaction energy is 

E(-£ n + £ n >) = - Jcos(2tt(-4 + l n ,)/N) (2) 

the site variable is £ n = 0, 1, ...N — 1, and (nn') denotes nearest-neighbor pairs. Equivalently, 
E = —JS n ■ S n >, where the angles of the (classical) spins take on the discrete values S n = 
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(cos# n , sm9 n ) with 9 n = 2-K£ n /N . We use the notation K = (3 J and u = e~ K ^ 2 . We consider 
the case N = 6 here because it provides a simple exactly solvable example of an intersection 
point with odd n c on a complex-temperature phase boundary B. The (reduced) free energy 
is f — —f3F = limjv s ->oo N' 1 hiZ in the thermodynamic limit. 

For d = 1 dimension, one can solve this model exactly, e.g., by transfer matrix methods. 
One has 
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Z = Tr(T N *) = J>f' (3) 

where the Xj, j = 1, ...6 denote the eigenvalues of the transfer matrix T defined by T nn i = 
(£ n \exp(—/3E(£ n ,£ n i))\£ n t), and we use periodic boundary conditions. It is convenient to 
analyze the phase diagram in the u plane. For physical temperature, phase transitions are 
associated with degeneracy of leading eigenvalues [16|]. There is an obvious generalization 
of this to the case of complex temperature: in a given region of u, the eigenvalue of T 
which has maximal magnitude, X max , gives the dominant contribution to Z and hence, in 
the thermodynamic limit, / receives a contribution only from X max : / = ln(X max ). For 
complex K, f is, in general, also complex. The CT phase boundaries are determined by the 
degeneracy, in magnitude, of leading eigenvalues of T. As one moves from a region with 
one dominant (i.e., leading) eigenvalue X max to a region in which a different eigenvalue X' max 
dominates, there is a non-analyticity in / as it switches from / = \ia(X max ) to / = ln(X' max ). 
The boundaries of these regions are defined by the degeneracy condition among dominant 
eigenvalues, \X max \ = \X' max \. These form curves in the u plane. Note that the free energy 
and the conditions for CT phase boundaries are the same if some eigenvalues occur with 
finite multiplicity greater than one. To see this, assume Xj occurs rij times, where rij is 
finite. If Xj is nonleading, the result is obvious; if Xj is leading, the result follows because 
limjv^oo iV s -1 hi(rijXj) = In Xj for any finite (nonzero) rij. This is relevant here because two 
of the eigenvalues of T occur with multiplicity 2. 

Of course, a ID spin model with finite-range interactions has no non-analyticities for any 
(finite) value of K, so that, in particular, the ID Zq model is analytic along the positive 
real u axis. Because of the invariance of Z on a bipartite lattice under the transformation 
S e — > S e , S — > —S , J — ► — J, where e and o denote sites on even and odd sublattices, it 
follows that 

u — > l/u B invariant (4) 

(For a finite ID lattice with periodic boundary conditions, we preserve this invariance by 
using even N s .) Further, since the Xj are analytic functions of u, whence Xj(u*) = Xj(u s )*, 
it follows that the solutions to the degeneracy equations defining the boundaries between 
different phases, |A»| = \Xj\, are invariant under u — > u* . Hence, B is invariant under u — > u* . 
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Although the model has a physical phase structure consisting only of the -^-symmetric, 
disordered phase, its complex-temperature phase diagram is nontrivial and exhibits a number 
of interesting features. An important property of this model is that the interaction energies 
E(A£), where A£ = i n — £ n i, are all rational multiples of each other: E(A)/J = —1, —1/2, 
1/2, and 1 for M = 0, ±1, ±2, and ±3, respectively, with E(6 - Al) = E(A£). Hence, all 



of the Boltzmann weights are powers of u 
has the Toeplitz form 

/ 1 
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and we set A34 = A3 = A4 and A56 = A5 = Xq. 

We find the complex-temperature phase diagram shown in Fig. 1, consisting of five 
different phases. In the caption, each phase is labelled by the eigenvalue (or identically equal 
eigenvalues) which is (are) dominant within it. The first phase is the complex-temperature 
extension (CTE) of the physical Z 6 -symmetric, paramagnetic (PM) phase, denoted R2 since 
A2 is dominant in this phase. The other four regions have no overlap with any physical phase 
and are thus of O ("other") type in our previous notation fll| . These include two phases -R34 



and -R34 in which A34 is dominant; these are related to each other by complex conjugation 
and occupy regions shaped like truncated lobes in the Re(u) > half-plane. Next, there is 
the phase R\ which occupies most of the left-hand plane Re(u) < together with a narrow 
portion to the right of the vertical axis Re(u) = 0. The fifth phase, R^ comprises an oval-like 
region centered on the negative real u axis in which A56 is dominant. A striking feature of 
the CT phase diagram is the presence of four intersection points on the CT boundary B, 
at u ~ 0.11 + 0.512 (given more precisely below), its inverse, and their complex conjugates. 
This is quite different from the CT phase diagram for the ID higher-spin Ising model [ B? |, 
where the curves (and, for half-integral spin, the semi-infinite line segment on the negative 
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Figure 1: Phase diagram of the ID Z§ clock model in the complex u plane. Phases are 
(a) i? 2 m the Re{u) > half-plane, extending in from the far right; (b) i? 34 occupying the 
truncated lobe-like region in the northeast quadrant, and its complex conjugate phase, R^; 
(c) Ri occupying most of the Re(u) < region, and (d) the oval-like phase R^ centered on 
the negative real u axis. 



real u s axis) only meet at the origin of the complex u s plane, i.e., at K = oo, where the 
model has an FM critical point and a formal essential singularity. This intersection point 
has n c = 4s 2 . (Given the u — > 1/win the latter model, in the u~ x plane, the curves also meet 
at the AFM critical point at uj 1 = 0.) A further difference is the presence in the Z G model 
of CT phases of finite extent in the u plane, viz., R34, R% 4 , and R§§. In the ID spin s Ising 
model, all CT phases extend from the wedge which they occupy in the vicinity of the origin 
outward to the circle at infinity in the u s plane. A third comparison may be made with the 
CT phase diagram for the ID g-state Potts model, which has no intersection points on B and, 
for q > 3, just two phases in the e~ K plane, the CT extension of the PM phase, of infinite 
extent, and an O phase of finite extent, bounded by the circle e~ K = (—1 + e llJJ )/(q — 2) [ |TU| . 
(The q = 2 Potts model is equivalent to the spin 1/2 Ising model.) 

To discuss the CT phase diagram of the ID Z 6 model, we consider first the oval boundary 
of the phase R^. This is the simplest portion of B because it does not contain any intersection 
points. It is the locus of solutions to the degeneracy condition |Ai| = [ A56 1 where these 



alternate as the dominant eigenvalues. Let u = re 10 . Then 
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\56| 



r + 2(1 + r 2 ) cos 9 + 4r cos 29 



(10) 



Note that in obtaining the equation in r and 9 above from the condition |Ai| 
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have removed the common factor 11 
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2 , which is irrelevant for £>, since at u = 1 the 
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condition that Ai and A56 are dominant is not satisfied (indeed, both vanish there). Further, 
in eq. fllCf) , the solution r = 0, i.e. u = 0, is not relevant for B, since A56 is not a dominant 
eigenvalue at this point. Equating the remaining factor to zero in the region where Ai and 
A56 alternate as dominant eigenvalues, one obtains the oval curve shown in Fig. 1. (If one 
relaxes the condition that Ai and A56 be dominant eigenvalues, then one obtains also a curve 
in the Re(u) > half-plane passing through u = 0, extending in "northeast" and "southeast" 
directions; however, this is not relevant for B.) Equating the factor in square brackets in eq. 
(pi]) to zero and solving, one gets 

1 



COS#i- 



56,± 



( r -l 



-) ± Vr- 2 + 26 + r 2 



(11) 



Observe the r — ■> 1/r symmetry in (|ll]), in accord with (|4]). The — sign gives the oval 
while the + sign gives the above-mentioned curve in the Re(u) > region, which does not 
contribute to B. Setting 9 = tt in (|T0|) yields the solutions r = 1/2 and r = 2, i.e. u = — 1/2 
and u = —2; these are the points where the oval boundary crosses the negative real u axis. 
Setting r = 1 in the — case of eq. (|ll|) yields the values of 9 at which the oval boundary 
crosses the unit circle \u\ = 1, which are given by ±0 O , where 

9 a = arccos^- 1 = 155.705° (12) 

The corresponding points are 

1 r -(l + V7) ± (8 - 2v / 7) 1/2 l ~ -0.911 ± 0.411z (13) 



As one travels along the boundary curve of the oval, 9 starts at the value tt at u = —1/2, 
decreases to the value 9 Q as u crosses the unit circle, and then increases back to 9 = tt at 
u = —2, and similarly with the complex conjugate points. 

The portion of B in the Re(u) > half-plane is more complicated and consists of several 
parts, each of which represents a degeneracy condition of the form |Aj| = |Aj| where Aj and 
Xj alternate as dominant eigenvalues. The part of B near the origin (and hence also, given 
the symmetry (f|), the part farthest out from the origin) is a curve consisting of the solution 
to the condition 

|Ai| = |A 2 | =^ r[2(l + r 6 )cos# + 4r 3 cos2# + r 3 cos4#] = (14) 

where Ai and A2 are dominant eigenvalues. For r ^ 0, we analyze the second factor in square 
brackets. For small r this reduces to the condition that cos# = 0, i.e., 9 = ±7r/2, so that this 
curve passes through the origin in a vertical direction. (Equation (0) also has a solution 
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consisting of a small oval which crosses the negative real u axis at — 2 ±1//3 , but this is not 
relevant to B since Ai and A2 do not alternate as dominant eigenvalues in this region.) 

The left-hand boundaries of the lobe regions -R34 and R\ A are given by the solution to 
the condition 

I Ai| = IA34I r[r - 2(1 ±r 2 )cos#] = (15) 

i.e., cos 9 = r/[2(l + r 2 )], where these eigenvalues are dominant. The corresponding curve 
crosses the unit circle at 9 = arccos(l/4) ~ ±75.52°. 

The right-hand boundaries of these lobe regions -R34 and are given by the solution 
to the condition 



I A* I = IA 



31 



3r(l + r 4 ) + 2(l + r 6 )cos# + 10r 3 cos2#±6r 2 (l±r 2 ) cos3#±4r 3 cos 49 



= 
(16) 

where A2 and A34 alternate as dominant eigenvalues. This curve crosses the unit circle at 
cos# = a 1 / 3 + (3/8)a- 1 /3 _ y 2 where a = (9 + 3V3)/32, i.e., 9 = ±41.03°. 

As noted above, there are four intersection points, each with n c = 3, on B. We denote the 
one closest to the origin in the northeast quadrant as u t = r t e t6t ; the others are then l/u t , 
and their complex conjugates, u* t and l/iz£. The point u t occurs where the three eigenvalues 
Ai, A2, and A34 are all degenerate and dominant, i.e., where the phase boundaries defined by 
the conditions |Ai| = |A 2 |, |Ai| = |A 34 |, and \X 2 \ — |A 34 | all meet. This point is given by the 
solution r t = 0.518684... to the equation 

2(r -6 ± r 6 ) — 18(r -2 ± r 2 ) — 31 = (17) 

This equation is manifestly symmetric under the symmetry (0), and has two positive real re- 
ciprocal roots, of which the smaller is the one which we have listed above. The corresponding 
angle is 9 t = 78.207851°, so that 

u t ~ 0.1059993 + 0.5077375i (18) 

A general characteristic of ID spin models with short-range forces is that the points 
K = ±00 are critical points, so that in the present case, the phase boundary B passes 
through the respective origins in the u and u~ l planes. This is reminiscent of how B passed 
through the origin of the u s plane for the ID spin s Ising model and the e~ K plane for the 
ID g-state Potts model 0. 

In order to gain further insight into the components of B, one may determine the full 
locus of points comprising the solution of the degeneracy condition |A»| = \Xj\ without the 
constraint that A, and Xj be dominant eigenvalues. One can then see how subsets of this 
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Re(u) 

Figure 2: Solutions with Re(u) > of the conditions (a) |Ai| = | A2 1 (solid curve); (b) 
I Ai| = | A34I (dashed curve); and (c) |A 2 | = |A 34 | (dot-dashed curve). 

locus form the various components of B. In Fig. 2 we show the portions of the solutions to 
the |Aj| = \Xj\ equations for = (1,2), (1,34), and (2,34) which are involved with the 
finite- if intersection points. The loci of solutions to | A2 1 = | A34 1 and |Ai| = | A56 1 also contain 
curves which pass through the origins of the u and u^ 1 planes; we omit these from Fig. 2 
since they do not contribute to B except at the respective origins of the u and u~ l planes, 
which are the zero-temperature FM and AFM critical points of the model. The intersection 
point u t , together with its inverse and their complex conjugates, are seen to be the (finite- If) 
points at which all of the curves |A;| = \Xj\ for = (1,2), (1,34), and (2,34) cross. While 
a solution curve of |Aj| = \Xj\ for a given (i, j) continues smoothly through u t , the actual CT 
phase boundary B has a discontinuous tangent there and heads off in a different direction 
as one moves away from u t , since it is then determined by a different degeneracy condition 
|Aj| = |Afc|. 

We generalize our results as follows: for a ID discrete spin model with short-ranged 
interactions, so that the transfer matrix has a finite number of eigenvalues, the CT phase 
boundary B has an intersection point with n c curves meeting if and only if n c (in general 
distinct) dominant eigenvalues become degenerate in magnitude at this point. 

Of course, there are differences between d = 1 and d > 2 models. For d > 2, T has an 
infinite number of eigenvalues, rather than the finite number which it has for d — 1. Further, 
since d > 1 is above the lower critical dimensionality de. c .d. = 1 of a usual discrete spin model, 
two eigenvalues which are distinct in the symmetric phase may be degenerate in the phase 
with spontaneously broken symmetry and ferromagnetic (FM) or antiferromagnetic (AFM) 
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long-range order, again in contrast to the situation in d = 1, where there is no symmetry- 
breaking phase for finite K. However, these differences do not prevent the ID Z$ model from 
having the same feature, viz., an intersection point on B with higher-dimensional 
models such as the spin 1/2 3D Ising model. 

Indeed, we have previouly pointed out some intriguing connections between certain fea- 
tures of the respective CT phase boundaries of the higher-spin Ising model in ID and 2D 



fl8| , |19] . These included the observation that (exactly known) angles at which phase bound- 
aries crossed the unit circle in the ID case were consistent with being equal to values of 
analogous angles in the 2D model. Remarkably, we find the same kind of connection here. 
The physical phase structure of the 2D Z 6 model consists of (i) a disordered, Z 6 -symmetric 
high-temperature phase, (ii) a Z 6 -symmetric intermediate-temperature phase with algebraic 
decay of correlations, and (iii) a broken-symmetry low-temperature phase with FM (AFM) 
long-range order for J > (J < 0) |j20| . From the CT zeros of the partition function for 



the Zq model calculated on a 6 x 7 lattice [21|, which we have confirmed in an independent 



calculation, one infers that in the "northwest" and "southwest" regions of the u plane, in 
the thermodynamic limit, B has intersection points with n c = 4, n t = 2, and m = 2 at 
e ±i9 °, where 9 is the same angle as we have calculated in eq. ([12]) . This suggests that one 
can gain valuable information about certain features of the complex-temperature properties 
of a higher-dimensional spin model from the CT properties of the corresponding ID model, 
reminiscent of the fact that e expansions starting from d = G^. c .d. gave similar insight into 



the physical phase transition in models for d > dg, c ,d. [22 
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